Subgap anomaly and above-energy-gap structure in chains of diffusive SNS junctions 
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We present the results of low-temperature transport measurements on chains of superconductor- 
normal-constriction-superconductor (SNS) junctions fabricated on the basis of superconducting PtSi 
film. A comparative study of the properties of the chains, consisting of 3 and 20 SNS junctions in 
series, and single SNS junctions reveals essential distinctions in the behavior of the current-voltage 
characteristics of the systems: (i) the gradual decrease of the effective suppression voltage for the 
excess conductivity observed at zero bias as the quantity of the SNS junctions increases, (ii) a rich 
fine structure on the dependences dV/ dl — V at dc bias voltages higher than the superconducting gap 
and corresponding to some multiples of 2A/e. A model to explain this above- energy- gap structure 
based on energy relaxation of electron via Cooper-pair-breaking in superconducting island connecting 
normal metal electrods is proposed. 
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In the past few years, the mesoscopic systems, consist- 
ing of a normal metal (N) or heavily doped semiconduc- 
tor being in contact with a superconductor (S), have at- 
tracted an increased interest mainly because of the rich- 
ness of the involved quantum effects |J. The key mech- 
anism governing the carrier transport through the NS 
contact is the Andreev reflection In this process, an 
electron-like excitation with an energy e smaller than the 
superconducting gap A moving from the normal metal to 
the NS interface is retro-reflected as a hole-like excitation, 
while a Cooper pair is transmitted into the superconduc- 
tor. This phenomenon is the basis of the proximity effect, 
which generally implies the influence of a superconductor 
on the properties of a normal metal when being in electri- 
cal contact. The consequences of Andreev reflection on 
the current- voltage characteristics (CVC) of a NS junc- 
tion are studied in detail in the so-called BTK theory || , 
which describes the subgap current transport in terms of 
ballistic propagation of quasiclassical electrons through 
the normal-metal region, accompanied by Andreev and 
normal reflections from NS interface. The probability 
of Andreev reflection and normal reflection are energy- 
dependent quantities and the relation between them de- 
pends on the barrier strength at NS interface, that is 
characterized by a parameter Z ranging from for a per- 
fect metallic contact to oo for a low transparency tunnel 
barrier. For a perfect contact (Z = 0), probability of 
Andreev reflection is equal to unity for particles with 
an energy e smaller than the superconducting gap A and 
the subgap conductance are found to be twice the normal 
state conductance thus demonstrating the double charge 
transfer. In the other limit, when Z — > oo, the conduc- 
tivity is very small. 

When a normal metal is placed between two super- 
conducting electrods another mechanism is involved in 
charge transfer. It is multiple Andreev reflection process 
(MAR). The concept of MAR was first introduced by 
Klapwijk, Blonder and Tinkham in order to explain 
the subharmonic energy gap structure observed as dips 



in the differential resistance dV/dl of SNS junctions at 
the voltages V n = 2A/ne (n = 1,2,3, •••). In all sys- 
tems, diffusive as well as ballistic, the MAR mechanism 
relies on the quasiparticles to be able to add up energy 
gained from multiple passages. This energy gain results 
in strong quasiparticle nonequilibrium. Calculations of 
the CVC, within the approach developed in Ref. ||], and 
taking into account MAR process was performed in || 
(OTBK theory). It was shown that MAR process results 
not only in SGS at voltages V n = 2A/ne but affects the 
general form of CVC at all voltages as well. Although 
BTK theory is not intended to describe diffusion trans- 
port in normal- metal region, a suitable value Z ~ 0.55 
of the barrier strenght gives result very close to those 
obtained by Green's-function method for microconstric- 
tion in the dirty limit Experimentally, in some cases 
BTK theory has been successful in providing a method 
to extract the junction reflection coefficients based on 
experimental CVC g§. 

At present, a large variety of SNS junctions (most of 
them are diffusive) is being fabricated and studied [0-0 . 
The investigations of these junctions are primarily fo- 
cused on the nonlinear behavior of the current-voltage 
characteristics, which exhibit the anomalous resistance 
dip at zero bias [zero bias anomaly (ZBA)], the subhar- 
monic energy gap structure (SGS), etc. These experi- 
ments have revealed a number of peculiar properties of 
diffusive SNS junctions, the explanation of those being 
beyond the ballistic BTK and OTBK theories. Nowa- 
days the properties of diffusive SNS junctions are con- 
sidered to be determined not only by the parameters of 
the NS interface, but nonlocal coherent effects as well, 
namely (i) the superposition of multiple coherent scat- 
tering at the NS interfaces in the presence of disorder 
(so-called reflectionless tunneling [p^[ ) and (ii) electron- 
electron interaction in the normal part. The latter is 
one of the important points of a recently developed "cir- 
cuit theory" when applied to diffusive superconducting 
hybrid systems [pi. Within this approach, which is 
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based on the use of nonequilibrium Green's functions, the 
electron-electron interaction induces a weak pair poten- 
tial in the normal metal. It results not only in a change 
of the resistance, but also in a non-trivial distribution of 
the electrostatic and chemical potentials in the structure, 
which implies non-local resistivity. Following the spirit 
of Nazarov's circuit theory Bezuglyi et al. [|l7| have de- 
veloped a consistent theory of carrier transport in long 
diffusive SNS junctions with arbitrary transparency of 
NS interfaces. Although much work both theoretical and 
experimental has been done on single SNS junctions, it is 
a challenge to fabricate and carry out comparative mea- 
surements on multiply connected SNS systems. In this 
paper, we present the results of low-temperature trans- 
port measurements on chains of SNS junctions and on 
single SNS junctions and perform the comparative anal- 
ysis of their properties. 

The design of our samples is based on the fabrication 
technique of SNS junctions that has been recently pro- 
posed and realized by us for preparation of single SNS 
junctions and two-dimensional arrays of SNS junctions 
[l4[ . The main idea employed in those experiments 
was to use superconductive and normal region made of 
the same material. The point is the suppression of su- 
perconductivity in the submicron constrictions made in 
an ultrathin polycrystalline PtSi superconducting film 
by means of electron beam lithography and subsequent 
plasma etching. Here we use the same fabrication tech- 
nique to prepare chains of SNS junctions. 

The original PtSi film (thickness - 6 nm) was formed 
on a Si substrate. The film was characterized using Hall 
bridges 50 /im wide and 100 /im long. The film had 
a critical temperature T c = 0.56 K. The resistance per 
square was 104 Q. The carrier density obtained from 
Hall measurements was 7 • 10 22 cm -3 , corresponding to 
a mean- free path I = 1.2 nm and a diffusion constant 
D = 6 cm 2 /s, estimated using the simple free-electron 
model. 

Single SNS junction is a constriction between two holes 
made in the film and placed in the Hall bridge. The 
hole diameters are 1.7 /im and the distance between cen- 
tres of holes is 2.1 /im, resulting in a width of narrowest 
part of the constriction of 0.4 /im. A scanning electron 
micrograph and a schematic view of a sample are pre- 
sented in Fig. la, b. The chains of SNS junctions are a se- 
ries of constrictions connected by the islands of the film 
(Fig.lc,d). They are designed to provide the possibility 
of the comparative study. The dimensions of constric- 
tions are identical with those of the single SNS junction, 
and the characteristic size of the islands of the film is 
~ 1.3 /im. As the constrictions are not superconducting 
we have a chain of SNS junctions (Fig.le). 

The measurements were carried out with the use of 
a phase sensitive detection technique at a frequency of 
10 Hz that allowed us to measure the differential resis- 
tance (dV/dl) as a function of the dc current (I). The 



ac current was equal to 1-10 nA. Figure 2 shows typi- 
cal dependences of dV/dl-V for the samples with single 
constriction (1D-1) and chains, consisting of 3 (sample 
1D-3) and 20 (sample 1D-20) SNS junctions in series. 
On the abscissa the value of voltage obtained by numeri- 
cal integration of the experimental dependences dV/ dl-I 
and afterward divided by the number of SNS junctions is 
plotted. Measured differential resistance is also divided 
by the number of SNS junctions (by 3 for structure 1D- 
3 and by 20 for structure 1D-20). Thus for the chains 
data presented in Fig. 2 are dependences the differential 
resistance as a function of dc bias voltage per SNS junc- 
tion in the average. As is seen at high voltage for chains 
the average resistances of SNS unit of chains obtained 
in such a way are close to each other and to the value 
dV/dl of the structure with single SNS junction. A com- 
mon feature for all structures studied is a minimum of 
the differential resistance at zero bias voltage. It points 
out the high transparency of NS interfaces in these sam- 
ples. The data for single SNS junction (1D-1) exhibit 
a behavior very similar to that reported in our previous 
work [[l3| . In Ref. |l7[ the differential resistance of long 
diffusive SNS junction with perfect NS interfaces at zero 
voltage was estimated by the following expression: 

dV/dI(0) = R N (1 - 2.64^/L), (1) 

where £,n = ^/TiD/2nkT is the decay lenght for the pair 
amplitude, L is the lenght of the normal-metal region, 
and Rn is its resistance at T > T c . The expression 
(1) has obtained on condition that the inelastic lenght 
l e exceeds the junction lenght L. The energy relaxation 
is described by the time r e _1 = ire 2 /(ShEp). This is 
the time it takes for a "hot" quasiparticle with energy e 
much larger than temperature T to thermalize with all 
the other electrons. For the excitation energy e = 2A 
(A/e - 270 /iV) we obtain l e = ^/D7 e ~ 4.6 /im. The 
magnitude of l e really exceeds the characteristic lenght of 
the constriction. As Rn we take the difference between 
the resistance of the whole structure with constriction 
and the resistance of the original film without constric- 
tion measured between the same probes at T > T c . It 
gives Rn ~ 530 Q. Assuming L ~ 1 /im, we find from 
(1) at T = 100 mK the differential resistance at zero 
bias dV/dI(0) ~ 0.SR N = 420 n. As is seen in Fig.2, 
this value is close fit to measured dV/dI(0) for the sam- 
ple with single constriction 1D-1. A clear-cut distinction 
between dependences dV/dl-V of the single SNS junc- 
tion and the chains consists in non-monotone behavior 
of dV/dl-V characteristics of the latter. The differen- 
tial resistance of the chains has a minimum at zero bias 
voltage and shows a maximum at a finite bias voltage of 
about 250 /iV for the sample 1D-3 and 30 /iV for 1D- 
20. The second feature is the gradual decrease of the 
effective suppression voltage for the excess conductivity 
observed at zero bias as the quantity of the SNS junctions 
increases. The same behavior - non-monotonic dV/dl-V 
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characteristics and the decrease of the effective suppres- 
sion voltage for the excess conductivity in comparison 
of that in single SNS junctions - has been observed in 
two-dimensional arrays of SNS junctions jL4|. 

The most unexpected result obtained for chains of SNS 
junction is the presence of a fine and fully symmetri- 
cal structure as dips in the differential resistance at high 
biases (Fig. 3a, b). The positions of these dips approxi- 
mately correspond to 4, 6, and 12 multiples of the su- 
perconducting gap. As is seen in Fig. 3c, the temperature 
dependence of the positions of these dips actually reflects 
the dependence A(T). To our knowledge the results pre- 
sented in this paper is the first investigation of chains of 
SNS junctions and the first observation such above en- 
ergy gap structure (AGS). To offer an explanation of the 
AGS let us consider a current driven in the normal re- 
gion at high voltages. As the inelastic lenght more than 
the normal metal region, at finite voltage quasiparticle 
passing through Ni between two superconductors Si and 
S2 (Fig. 4b, d) results in a nonequilibrium energy distribu- 
tion of quasiparticles. For a start we address this problem 
to the ballistic OTBK theory |^|. In this approach the 
quasiparticles are divided into two subpopulation which 
depend on their direction of motion /_>(22) and /<_(£"), 
with all energies measured with respect to the local chem- 
ical potential. The current through the junction is 
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In Fig.4a,c we have plotted /-(£), /<-(#), and [/-(£)- 
/<_(£")], assuming a parameter Z = 0.55. The volt- 
age across the junction equals 4A/e and 6A/e in Fig. 4a 
and Fig. 4c, respectively. The function [/_>(£■) — /<_(£")] 
is depicted by the blackness (high values are black) in 
Fig. 4b, d, with the energy scale being plotted vertically. 
These calculations clearly show that the nonequilibrium 
distributions are sharply peacked at the four gap edges. 
What happens to quasiparticle injected to the supercon- 
ducting electrode S2 and having the energy more than 
2 A measured with respect to the gap edge of S2? They 
can decay under spontaneous phonon emission into states 
of lower energy [ 18[ . Phonons emitted in the process 
of energy relaxation of injected quasiparticles can be re- 
absorbed via Cooper-pair-breaking. Depending on the 
volume of finite states, the probability of this process is 
non-monotone function of the primary energy of injected 
quasiparticles. It takes the peak values at the energies 
multipled to 2 A, for density of quasiparticle states to 
have the singularity at the gap edge. As a result one 
quasiparticle injected to S2 with energy high the gap edge 
by m-2A can cause up to l+2m quasiparticles injected to 
N2 (Fig. 4b, d) and consequently the increase of the differ- 
ential resistance. The realization of this scenario requires 
the lenght of energy relaxation to be comparable to di- 
mensions of superconductor. It is really fulfilled in the 



SNS systems studied. Proposed mechanism supposes the 
sistem consisting of two SNS junctions in series to be 
sufficient in order to observe the AGS. 

In summary, we performed the first comparative study 
of the properties of the chains and single SNS junctions. 
Our experiments reveal essential distinctions in the be- 
havior of the current- voltage characteristics of these sys- 
tems. A detailed quantitative analysis needs to take into 
account the contributions of nonlocal coherent transport 
and an effect of quasiparticle injection in both the su- 
perconducting and normal regions of chains of the SNS 
junctions. 
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FIG. 1. (a) Scanning electron micrograph of the sample 
with a single constriction formed by electron beam lithogra- 
phy and subsequent plasma etching of the 6 nm PtSi film 
grown on a Si substrate, (b) Schematic view of a junction 
(not to scale) showing the layout of the constriction in the 
Hall bridge, (c) Scanning electron micrograph of the sample 
with 20 constrictions, (d) SEM subimage of the sample rep- 
resented in (c). (e) The layout of a chain of SNS junctions 
showing the dimensions of the structure. Regions of the nor- 
mal metal constrictions are dark, and the superconducting 
islands are light gray. 
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FIG. 2. Differential resistance vs dc bias voltage for sam- 
ples: (1D-1) with single constriction, (1D-3) with three con- 
strictions in series, and (1D-20) with twenty constrictions in 
series. T — 100 mK. Differential resistance and dc bias volt- 
age are divided by 3 and 20 for samples 1D-3 and 1D-20, 
respectively. 
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FIG. 3. (a) Temperature evolution of the differential resis- 
tance of the sample 1D-20 as a function of the average bias 
voltage falling at one SNS junction. All traces except the 
lowest trace are shifted up for clarity, (b) The same for the 
sample 1D-3. The differential resistance reveals symmetrical 
minima at voltages exceeding 2 A. The arrows indicate the 
above- energy- gap structure (AGS), corresponding to integer 
multiples of the 2 A. The temperature dependence of the AGS 
at eV = m • 2 A divided by m for both samples is depicted 
in (c) by symbols, and compared to the BCS temperature 
dependence of the gap (solid line). 
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FIG. 4. Nonequilibrium distributions in normal region, for 
SNS junction with Z = 0: (a,b) at eV = 4A; (c,d) at eV = 6A 
and T = 0. In (b) and (d) the function [/-(£) - /<-(£)] is 
depicted by the blackness (high values are black) . The energy 
scale is plotted vertically. 
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